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We determine the concentration-mass relation of 19 X-ray selected galaxy clusters from the CLASH survey in 
theories of gravity that directly modify the lensing potential. We model the clusters as NFW haloes and fit their 
lensing signal, in the Cubic Galileon and Nonlocal gravity models, to the lensing convergence profiles of the 
clusters. We discuss a number of important issues that need to be taken into account, associated with the use of 
nonparametric and parametric lensing methods, as well as assumptions about the background cosmology. Our 
results show that the concentration and mass estimates in the modified gravity models are, within the errorbars, 
the same as in ACDM. This result demonstrates that, for the Nonlocal model, the modifications to gravity are 
too weak at the cluster redshifts, and for the Galileon model, the screening mechanism is very efficient inside the 
cluster radius. However, at distances ~ [2 — 20] Mpc /h from the cluster center, we find that the surrounding 
force profiles are enhanced by ~ 20 — 40% in the Cubic Galileon model. This has an impact on dynamical mass 
estimates, which means that tests of gravity based on comparisons between lensing and dynamical masses can 
also be applied to the Cubic Galileon model. 


I. INTRODUCTION 

Over the past decade, the ACDM paradigm has established 
itself as the standard model of cosmology. Most of the mat¬ 
ter in this model is in the form of cold dark matter (CDM), 
and a cosmological constant. A, plays the role of the dark en¬ 
ergy that is responsible for the late-time accelerated expan¬ 
sion of the Universe. Photons, massive neutrinos and baryons 
make up the rest of the energy budget. The final ingredient is 
Einstein’s theory of General Relativity (GR), which describes 
the gravitational interaction between these components. Al¬ 
though it is true that this model is in remarkable agreement 
with most of the cosmological data gathered to date m, L 
is also the case that, in some aspects, the model still lacks 
compelling theoretical and observational support. Perhaps the 
most worrying of these shortcomings is the unnaturally small 
value of A compared to the predictions from quantum field 
theory. Another problem relates to the extrapolation of the 
regime of validity of GR from the Solar System (where it has 
been very well tested O) to cosmological scales, where there 
is still a shortage of stringent model-independent tests of grav¬ 
ity. These two problems of ACDM have been fuelling inter¬ 
est in cosmological models with modified theories of gravity, 
which is now a well developed branch of cosmology on both 
the theoretical ElSl and observational EE) levels. 

Here, we focus on observational determinations of galaxy 
cluster masses derived from lensing in theories of modified 
gravity. This is a topic that has not been extensively investi¬ 
gated in the literature. The reason for this, we believe, is his¬ 
torical as many of the first modified gravity models to be com¬ 
pared to observations were models like f(R) liTOl or Dvali- 
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Gabadadze-Porrati (DGP) [QT| gravity, which do not mod¬ 
ify the lensing potential directly through a modified Poisson 
equation. In these models, which are conformally equivalent 
to scalar-tensor theories, lensing mass estimates are automat¬ 
ically the same as in GR, whereas dynamical mass estimates 
are not mm. This led to the development of a number of 
tests of gravity. For instance. Refs. mm compared galaxy 
lensing masses from strong lensing with dynamical masses 
inferred from stellar velocities to probe gravity on kpc scales. 
Also, Refs. liT6Hl8l developed methods to probe the modified 
dynamical potential in f(R) and DGP models in the infall re¬ 
gions around massive clusters, given the lensing mass. More 
recently. Ref. m used comparisons between the X-ray sur¬ 
face brightness and lensing profiles of galaxy clusters to con¬ 
strain models like f(R) (see also Ref. lf20l ). The unmodi¬ 
fied lensing potential in these theories also allowed for clus¬ 
ter lensing masses to be used as a relatively model indepen¬ 
dent ingredient in observational tests of gravity. For instance, 
in the work of Refs. EfflZD, the authors used the fact that 
f{R) models modify the halo mass function to place observa¬ 
tional constraints using data from the abundance of clusters as 
a function of their lensing mass. 

Recently, there has been growing interest in models that 
also modify the way in which the lensing potential depends 
on matter density perturbations, such as Nonlocal gravity m~ 
|29) , Galileon gravity I30ll32l . massive gravity 0314391 . K- 
mouflage gravity {4014431 . and several other subclasses of 
Horndeski’s general theory fl4) . These modifications to the 
lensing signal give rise to a broader range of ways to test grav¬ 
ity. For example. Refs. Il45l l461 presented forecasts for future 
galaxy-galaxy lensing observations m showing characteris¬ 
tic signatures of some models of massive gravity. Models that 
change the lensing signal can also have a strong impact on 
the power spectrum of cosmic shear ll48l l49l and the lensing 
of cosmic microwave background (CMB) photons 18115011711. 
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Moreover, photons are a direct probe of the time evolution of 
gravitational potentials, which allows strong constraints to be 
placed upon modified gravity models via the integrated Sachs- 
Wolfe (ISW) effect 0 [52] E2. 

Another consequence of the modifications to the lensing 
potential is that this may introduce model-dependent system- 
atics in the estimation of cluster masses from lensing. The 
investigation of such biases and their connection with some of 
the above-mentioned tests of gravity is one of the main goals 
of this paper. We choose the Galileon and Nonlocal grav¬ 
ity cosmologies as working examples of models that directly 
modify the lensing potential. In the Galileon model, an ex¬ 
tra scalar degree of freedom gives rise to a fifth force at late 
times. The effects of this fifth force are appreciable on large 
cosmological scales, but are suppresed near massive bodies by 
means of a screening mechanism known as the Vainshtein ef¬ 
fect 85411568 . This model has been shown to provide a good fit 
to the CMB temperature, CMB lensing and baryonic acoustic 
oscillations (BAO) data in Refs. El El. Nonlinear structure 
formation in this model has been studied in Ref. 85 81 using the 
spherical collapse model and in Refs. 15911601 using N-body 
simulation. Reference ED studied the properties of dark mat¬ 
ter haloes, which were used to develop a halo model of struc¬ 
ture formation for Galileon gravity. This model is, however, 
under observational tension as it predicts a negative sign for 
the ISW effect 10, which is at odds with recent observations 
8628 . In the case of the Nonlocal gravity model, the modifica¬ 
tions to gravity on cluster scales can be fully parametrized by 
a time-dependent effective gravitational strength. This model 
has no screening mechanism, but Refs. l26l l63l have shown 
that, if the background evolution can be neglected locally, 
then the model becomes compatible with Solar System tests 
of gravity. Linear structure formation has been studied by 
Refs. I1271 [28l and Ref. 8291 performed the first N-body simu¬ 
lations of Nonlocal gravity cosmologies, which were used to 
study halo properties and to also construct a halo model for¬ 
malism. In addition to the modified gravitational law, both 
the Galileon and Nonlocal gravity models also modify the ex¬ 
pansion rate at late times. This is different from models such 
as f{R) or DGP gravity for which the expansion rate can be 
tuned to match that of ACDM. 

To estimate lensing masses in Galileon and Nonlocal cos¬ 
mologies, we model galaxy clusters as Navarro-Frenk-White 
(NFW) haloes (64), and fit the predicted lensing convergence 
signal to the data obtained from weak and strong lensing 
observations for 19 X-ray selected clusters from the Cluster 
Lensing and Supernova Survey with the Hubble Space Tele¬ 
scope (CLASH) 86514671 . Our analysis is similar to that per¬ 
formed in the context of GR in Ref. 1661 . In most of the 
current data analysis, one often makes model-dependent as¬ 
sumptions which may lead to results that are biased towards 
the assumed models. For example, the analysis of Ref. 8661 
assumes a fiducial ACDM background to compute angular 
diameter distances. Assumptions like these must be identified 
and carefully assessed before using the observations to test al¬ 
ternative models. On the other hand, the lensing data analysis 
of Ref. 8661 makes no assumptions about the mass distribu¬ 
tion of the clusters, which reduces the model dependency of 


the conclusions drawn from the observations, and makes it 
particularly well suited to tests of modified gravity. Given the 
subtle nature of some steps involved in the analysis of lensing 
data we shall pay special attention to them and explain how 
they can be taken into account. 

The rest of this paper is organized as follows. In Sec. [II] we 
describe the calculation of the lensing convergence of NFW 
haloes in ACDM, Galileon and Nonlocal gravity models, and 
in Sec.[III]we analyse the general predictions for each model. 
We describe our fitting methodology in Sec. IV where we 
comment also on the extra steps that one needs to take, in or¬ 
der to account for certain model-dependent assumptions made 
in the data analysis. In Sec. [V] we present our main results 
for the cluster lensing mass and concentration estimates in the 
three models we consider. We also discuss the link between 
the results found here and tests of gravity on large scales, 
namely those that were first designed for models that do not 


modify lensing. We summarize in Sec. VI 


II. LENSING EQUATIONS 

In this section we specify our notation and describe the cal¬ 
culation of the lensing quantities in the models we consider. 


A. Cluster lensing basics 

Throughout, we work under the commonly adopted setup 
for cluster lensing studies (see e.g. Refs. 86814708 for compre¬ 
hensive reviews). In particular, we consider a set of source 
galaxies at redshift z s , whose light gets deflected by a galaxy 
cluster at Zd- We use /’),/, D s and I),is to denote, respec¬ 
tively, the angular diameter distances between the observer 
and the lens, the observer and the sources, and the lens and 
the sources. We assume clusters are spherically symmetric 
and use the thin-lens approximation in which one neglects the 
thickness of the galaxy cluster compared to the much larger 
values of Dd, D s and D r / S . We also neglect the lensing distor¬ 
tions induced by foreground and background structures, com¬ 
pared to the lensing signal of the cluster. In our notation, 
r = sjx 2 + y 2 is a two-dimensional radial coordinate de¬ 
fined on the lens plane and with origin at the cluster center 
( x and y are cartesian coordinates); l denotes the optical axis 
(line-of-sight) direction, perpendicular to the lens plane, and 
with origin also at the cluster center; and R = \Jr 2 + l 2 , is a 
three-dimensional radial coordinate with origin at the cluster 
center. 

Light rays coming from the sources are deflected at the lens 
position by an angle a, which is related to the true (unob¬ 
served) angular position, /3, and the observed one, 9, by 


P = 9-a(9). (1) 

The local properties of the lensing signal are fully determined 
by spatial second derivatives of the scaled projected lensing 
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potential of the galaxy cluster, ip, which is given by 

n, 9 f Dda 

iP(9 = r/D d ) = * — $i en (r,Z)dZ, (2) 

D d D s c A J_ Dd 

where c is the speed of light and <l>| otl = ( ( 1 > + T) /2 is the 
total three-dimensional lensing potential. The two Newtonian 
potentials $ and T are defined by the perturbed Friedmann- 
Robertson-Walker (FRW) line element: 


ds 2 = (1 + 2'T/c 2 ) c 2 df 2 - a 2 (l - 2$/c 2 ) dx 2 , (3) 

where a = 1/(1 + z) is the cosmological scale factor (z is 
the redshift). The Jacobian matrix of the lensing mapping of 
Eq. |T} is given by 


= I” 1 ~ «-7i -72 

d0 [ [ -72 1 - « + 7iJ ’ 

where 

*(®) = l (ti$ m + d% y ) 1>, 

= ^V 2 t/>= 8 2 + 8 2 )ip 

is the lensing convergencJj/j and 

71 = \ (cl - dg y ) ip, 

72 = de x de y ip, 


(4) 

(5) 


( 6 ) 


are the two components of the complex lensing shear, |'y| = 
l7? + 7 |. The convergence is responsible for an isotropic 
focusing (or defocusing) of the light rays, whereas the shear 
field causes distortions in the shapes of the observed source 
galaxies. 

In lensing studies, one can split the analysis into the weak 
and strong lensing regimes. In the weak lensing regime, the 
directly observable quantity is the locally averaged complex 
ellipticity field in the lens plane, (e), which can be constructed 
from measurements of background galaxy shapes. At each 
point of the lens field, an average is taken over a number of 
nearby sources to smooth out the intrinsic ellipticity of the 
galaxies from that caused by the lens (see e.g. mmm/ 
Observationally, the field (e) is directly related to the reduced 
shear, g, (see e.g. Sec. 4 of Ref. E3) 


(e) 


9 = 7/(1 - «)• 


(7) 


The strong lensing regime takes place in the inner most re¬ 
gions of the lens. There, the lensing quantities k and 7 be¬ 
come large and the equations become highly nonlinear. As a 
consequence, highly distorted images like giant arcs or arclets 
and multiple images of the same background source can form. 
This happens close to the location of the critical curves of the 


1 The overbar on the V operator indicates that it is the two-dimensional 

Laplacian. Also, note that r = D d 8. 


lens, which are defined as the set of points on the lens plane 
where the lensing matrix, Eq. Q, becomes singular, i.e., 

det (df3/86) = (1 - nf - 7 2 = 0 . ( 8 ) 

Observationally, one identifies multiple images and giant-arcs 
to infer the position and shape of the critical lines. Then, 
given a theoretical prediction for n and 7 , one can check if 

det (8/3/89) vanishes at the location of the critical lines. 


B. Convergence in ACDM 

In GR, in the absence of anisotropic stress, $ = 'T, and as a 
result, the lensing potential is equal to the dynamical potential, 
4>i en = (<h + d/) /2 = $ = vk. Both satisfy the Poisson 
equation, 

V (r,i) $ (G 0 = 47rGp(r, l ), (9) 

where G is Newton’s constant and p(r, l) is the three- 
dimensional density distribution. The lensing convergence is 
obtained by integrating Eq. (j9j) along the line of sight. 


J ^ 7 , 1 ) 4 »(g l)dl = A-kG J p{r,l)dl = 47 rGS(r), ( 10 ) 

where E(r) is the surface mass density. The left-hand side of 
this equation can be manipulated as follows: 


J Vf r ^(r,l)dl = J V 2 <b(r, Z)dZ 

= V 2 | d>(r,l)dl + [dMr,l)}tD d « 


D s c 2 

DdsDd 


n{r = 9D d ), 


8?${r,l)dl 
V 2 r ip(r) 


«c 2 ^ 


where we have used Eqs. 0 and |[5} , and also the fact that 
first spatial derivatives of $ are negligible at cosmological 
distances away from the lens (thin-lens approximation). Com¬ 
bining Eqs. ( [TO} and ( [TT) yields 


4nGD ds D d ^_ £( 0 ) 

= -3 -A-’ 


( 12 ) 


where 


S c 


c 2 D s 
47rG D ds D d 


(13) 


is called the critical surface mass density for lensing. There¬ 
fore, in GR, the calculation of the lensing convergence re¬ 
duces to the evaluation of the projected two-dimensional den¬ 
sity profile of the cluster, which can often be done analytically 
(see Sec. HE below). The Hubble expansion rate in ACDM 
is given by 
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^ —^ — E 2 (z) — n m0 (l + 2 ) 3 + (1 — !l m o), (14) 

where Hq = lOO/i km/s/Mpc is the present-day Hubble ex¬ 
pansion rate (the subscript "0" denotes present-day values) 
and n m0 is the fractional background energy density of pres¬ 
sureless matter. Here, and throughout, we assume a spatially 
flat Universe and neglect the contribution to the expansion rate 
from radiation and massive neutrinos. Equation ( fl4| > is used 
in the calculation of the angular diameter distances that enter 
Eq. fB) and the relation between radial and angular scales, 
r = D d 9. 


C. Convergence in Galileon Gravity 

We focus on the cubic sector of the Galileon gravity model 
j30U32l . Its action is given by 


in which tp is the Galileon field, Ad 3 = MpiH 2 , Mpj = 
1/(87 tG) is the reduced Planck mass, □ = V M is the 
d’Alembert operator, C m is the matter Lagrangian density and 
Greek indices run over 0, 1, 2, 3. In flat spacetime, the above 
action is invariant under the Galilean shift d^ip d^tp + 
(where h t , is a constant four-vector). In spherical coordinates, 
the Poisson equation in the Galileon model leads to the follow¬ 
ing force law (see Refs. f58li598 for details about the deriva¬ 
tion) 


®:R _ GM(< R ) _ c 3 j. 2 dip,R 
R R 3 M 3 ^ R 


(17) 


where Sip(R) is the spatial perturbation of the Galileon 
scalar field about the background value, (p[z), M(< R) = 
47t f Q p(r')r , 2 dr' is the mass enclosed inside radius R and ,r 
denotes partial differentiation w.r.t. R. Equation differs 
from GR by having an extra source term which is governed by 


8 <P,r 

R 


4 Mpi 

3 fi 2 


(,') {< 


■*)■ 


+ 1-1 


GM{< R) 


R 3 


(18) 


S= d 4 xy/^g 


i^G~\ C2C2 ~\ C ^- U 


(15) 


where 1Z is the Ricci scalar, g is the determinant of the metric 
g jn/1 c 2 and c 3 are dimensionless constants, and C 2 and £ 3 are 
given by 


£2 = £3 = (16) 


with 


16 M, 


tV = 


pi 


9 f3if3 2 M 3 


GM(< r), 


(19) 


where B\ and fi 2 are two dimensionless functions of time. 
The quantity ry is a radial scale (often called the Vain¬ 
shtein radius) that roughly determines the distance from the 
halo/cluster center within which the modifications to gravity 


are suppressed. The combination of Eqs. (17 1 , (18 1 and its 
derivatives leads to 


d+R 

R 


1 1 _ 4 c 3 <p 2 

\ 3 MpiM 3 




GM{< R) 

W~ 


&,RR = G 


M{< R) lR 
R 2 


2M(< R) 
R 3 


3 c 3 /3i(p 2 

4 



3 ( r v /R ) 2 

2 y/(r v /R ) 3 + 1 


( 


r v \ 

rv ' R ~n) 


( 20 ) 

( 21 ) 


In the limit of large R, Eq. (20 1 can be written as (note that 


’v 


constant as R increases) 


with 


*,R _ y, M{< R) 
R hn R 3 


Glin = G 



2 c 3 ip 2 \ 

3MpiM 3 p 2 J ’ 


( 22 ) 


(23) 


being an effective "linearized" time dependent gravitational 
strength. On the other hand, when R becomes small, it is 
straightforward to show that (using, for instance, the NFW 
expressions shown below) 


„ GM{< R) 

R ~ R 3 ' ( ; 

That is, at sufficiently small radii, the force in the Galileon 
model is approximately the same as in GR, which is a direct 
consequence of the Vainshtein screening mechanism. 

In the Cubic Galileon model one also has that $ = 4/ l58l . 
which implies, like in GR in the absence of anisotropic stress, 
that the lensing potential is equal to the dynamical potential. 
We compute the convergence in the Cubic Galileon model 
by numerically integrating the three-dimensional Laplacian of 
the total potential, V 2< 1>, along the line of sight as 
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(n\ DdsDd 

Kie) = -D^f 

1 

“ 4ttGE, 


'-D d 




^ cTd ( $ -««( r >0 + 2 ^( r , 0 )dJ 


(25) 


(recall that in spherical coordinates, V 2 - —> ■. pR +2 '^). 

All that is left to specify is the time dependence of the 
background quantities that enter the above equations of the 
Galileon model. The time evolution of the Hubble parameter, 
(p, /3i and fa are given, respectively, by (81 


E{af = 

(p = 

Pi = 

p 2 = 


(2 mO a 3 + \J^mO a 6 +4(l^(2 m o) 


1 
2 

£H 0 /E(a), 
1 


6c 3 
2M 3 M pi 
V 2 


~ C 2 ~w ^ +2 m+^M* 


2c 2 


Pi- 


, (26) 
(27) 
(?8) 

(29) 




As in Ref. ®, we take C 2 = — 1 and the other two Galileon 
parameters are determined by fl m o as 


£ = V6(l-Ono), (30) 

C 3 = l/(6£). (31) 


D. Convergence in Nonlocal Gravity 

We take the model of Refs. (26l|27l as the representative 
case of a Nonlocal gravity model. The action is given by 


a -t» 


7^ - — 7^□- 2 7^ - C r 
6 


.(32) 


This (nonlocal) action can be cast in a more familiar (local) 
form given by (7314751 


A= 


16nG 


J dx 4 ^ 


1Z - -jrKS - £i (DU + 7^) 


-£ 2 (DS + C/)-£ m ], 


(33) 


where £ 3 and £ 2 are two Lagrange multipliers and we have 
introduced two auxiliary scalar fields, U = — CD 1 ^ and 
S = □ 2 77.. For completeness, we note that these two for¬ 
mulations are not equivalent and that care must be taken when 
matching the solutions associated with the above two actions 
(see e.g. Refs. (25l[75lj80l l for a discussion) 

Following Refs. (27] |29l, for the scales relevant for large 
scale structure formation, the modifed Poisson equation in this 
Nonlocal model can be written as 


(34) 


which takes the same form as in GR, Eq. (|9]), but with an ef¬ 
fective time dependent gravitational strength given by 


Glin = G 


1 - 




1 -1 


(35) 


where S is the background part of the field S. The time evo¬ 
lution of the background quantities in the Nonlocal model has 
to be obtained by numerically integrating the background dif¬ 
ferential equations (see e.g. Refs (27l[29l ). The parameter m 
in Eqs. (32] i and ( |33[ > controls the amount of dark energy in the 
Universe. In a flat Universe, the value of m is therefore deter¬ 
mined by the energy densities of the remaining matter species, 
which means this Nonlocal gravity model has the same num¬ 
ber of free parameters as ACDM. 

Just as in the cases of ACDM and Cubic Galileon gravity, 
in the Nonlocal model we also have that $ = d/ = <l>i eri , in 
the absence of anisotropic stress. Moreover, since the Pois¬ 
son equation in this model is obtained from GR by a simple 
rescaling of the gravitational strength, it follows that the con¬ 
vergence can also be computed analytically as 


k{9) 


G lin \ £(0) 
G ) £ c • 


(36) 


Note, however, that the Nonlocal expansion rate must be used 
in the calculation of the angular diameter distances that enter 

S c . 

Contrary to the case of Galileon gravity, the Nonlocal 
model does not possess a screening mechanism, which may 
raise some concerns about the ability of this model to pass 
Solar System constraints 0. For instance. Ref. (29l showed 
that the rate of change of the gravitational strength on cosmo¬ 
logical scales, Giin, if applied directly to Solar System tests, 
results in the model becoming inconsistent with current lunar 
laser ranging experiments ED. However, the time variation 
of Gun follows from the background expansion rate, and it is 
not clear what its impact is in the Solar System. In fact, the 
authors of Refs. I26ll63l have shown that if the cosmological 
expansion is neglected, i.e. if the spacetime about the Sun is 
perturbed Minkowskii (as opposed to FRW), then the model 
predictions become compatible with the current bounds. Here, 
we shall bear these discussions in mind, but focus instead on 
the model predictions for cluster scales, which are sufficiently 
large for one to consider the gravitational strength given by 
Eq. ((35]). 


E. NFW halo expressions 

In order to compute the lensing convergence in any of the 
cosmological models considered above, we need to specify 
the density profile of the lenses, which we model as dark mat¬ 
ter haloes with NFW density profiles (641 . 


Pnfw(-R, z) 


Ps{z) 

(R/r s )(l + R/r s ) 2 ' 


I* = 4:TrG e ff(z)p(r,l) 


(37) 
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This profile is fully specified by two parameters known as the 
scale radius, r s , and the characteristic density, p s . The mass 
enclosed inside radius R in a NFW halo is given by: 


M N fw(< R) = 47 rp s rf 


In (1 + R/r s ) 


R/r s 
1 + R/r s _ 

(38) 


We define halo masses in the usual way 



M a = — A Pc R 3 a , (39) 

where R A is the radius within which the mean density is equal 
to A times the critical density of the Universe at a given red- 
shift, p c = 3H 2 (z)/(8nG). Equating Ma = Mnfw(< Ra) 
one finds 


1 A 3 

Ps = ^t\p C C A 


In (1 + ca) - 


ca 


i -i 


1 + ca 


(40) 


where we have defined the concentration parameter 


CA = Ra/t s . (41) 

We take A = 200 (as it is standard in the literature since it 
is close to the overdensity at virial equilibrium) and instead of 
characterizing the NFW haloes by p s and r s , we use the equiv¬ 
alent and more common parametrization in terms of M 200 and 
C200- 

The surface mass density of a NFW halo admits an analyt¬ 
ical solution given by l68l [82l i83l 




DnfwO" = D d 6) = 


Pnfw(g l)dl = 


(42) 



six 1 — 1 


arctanh 


arctan 



x < 1 
x = 1 
x > 1 


where x = r/r s . The calculation of the lensing convergence 
in Galileon gravity also requires the evaluation of the gradient 
of Mnfw(< R), which is given by 

R 

AFnfw(< R),r = 47t p s r 3 s - - , 2 ■ (43) 

(r s + RY 


HI. MODEL PREDICTIONS FOR FIXED CLUSTER 
PARAMETERS 


To develop intuition about our results, we first analyse the 
model predictions for fixed cluster parameters, M^oo, C 200 , 
and cosmological matter density 1 \ nM . Figure [T] shows the 


FIG. 1. The upper and middle panels show, respectively, the time 
evolution of the Hubble expansion rate, H(z)/Ho, and of the ef¬ 
fective linear gravitational strength, Gu n , for the Galileon (red) and 
Nonlocal (blue) models, plotted as the relative difference to ACDM. 
The shaded band in these two panels indicates the redshift range 
spanned by the CLASH clusters analysed in this paper. The lower 
panel shows the density (solid) and enclosed mass (dashed) profiles 
for the Galileon (red) and Nonlocal (blue) models, plotted as the rel¬ 
ative difference to ACDM. The NFW parameters and redshift are 
M 200 = 5.0 x 10 14 Mq//i, C200 = 4.0, Zd = 0.5. To guide the eye, 
the shaded band in the lower panel indicates the radial scales outside 
R 200 in ACDM. 


time evolution of H(z)/H 0 (top panel) and G'| itl (middle 
panel) for the Galileon and Nonlocal models, plotted as the 
difference relative to ACDM. The lower panel of Fig.[l]shows 
the NFW density (solid) and enclosed mass (dashed) profiles 
for a halo at = 0.5, and with M 2 oo = 5.0 x 10 14 M 0 //i and 
C 200 = 4. For all models, the cosmological matter density is 
D m o = 0.27. We note, for completeness, that if h is absorbed 
into unit definitions (e.g. Mq//i for masses or Mpc//i for dis¬ 
tances), then our analysis becomes completely independent of 
its value, for all models. 
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FIG. 2. (Upper left) Radial profiles of the Laplacian (V 2< f>, solid) and gradient (4>, r , dashed) of the total gravitational potential, <f>, for the 
Galileon model (red), plotted as the relative difference to ACDM (black). The NFW parameters and redshift are M 200 = 5.0 x 10 14 Mq//i, 
C200 = 4.0, Zd = 0.5. Also shown are the predictions for a model called QCDM (green), which has the same background cosmology as 
the Galileon model, but with the force calculated as in GR. To guide the eye, the shaded band indicates the radial scales outside R 200 in 
ACDM. (Upper right) Same as the upper left panel but for the Nonlocal gravity model (blue). (Lower panels) Same as the upper panels but 
for the lensing convergence, k(0), and assuming z 3 = 1. In the lower panels we show d on the x-axis, which, in the different models, relates 
differently to physical distances due to the modifications to Dd- 


Figure |T| shows that the amplitude of the halo density pro¬ 
files near their center is lower in the Galileon and Nonlocal 
models, than it is in ACDM. This is because of the lower val¬ 
ues of p c {zd ) in these models. Specifically, in the inner-most 
part, the density profile of the halo becomes approximately 


Pnfw(R) 


p s r s 

R 


oc 


Pc(zd ) 2/3 

R 


(44) 


• —1/3 

as can be checked by noting that p s oc p c and r s oc p c , 
and recalling that we are assuming fixed M 200 and 0200 - As 
a result, if H(z) is smaller than in ACDM at Zd (upper pan¬ 
els), then so is pV^{zd) oc H A / 3 {zd), which effectively leads 
to a less dense halo. The same qualitative reasoning also ap¬ 
plies to the regime where the NFW density scales as oc R~ 2 . 
However, far from the halo centre, we have 


Pnfw( 7?) ~ oc (45) 

In this case, the dependence on p c {zd ) cancels out, and hence, 
all models have the same density values, as seen in the lower 
panel of Fig. |T](7?. > 10Mpc/A). The enclosed mass profiles 
show a qualitatively similar trend. Perhaps the only notewor¬ 
thy difference is that the M{< R) profiles for the different 


models do not agree at large radii. In this limit, the enclosed 
mass (Eq. ( |38| >) scales with radius as 

M(< R) oc [In (1 + R/r s ) - 1], (46) 

which retains a cosmological dependence (1 /r s oc p c ) in 
the logarithmic divergence of the mass. In particular, in the 
Galileon and Nonlocal models, p c is smaller than in ACDM, 
which implies that M(< R) is also smaller, as shown. 

Another important consequence of the modified Hubble ex¬ 
pansion rates in the Galileon and Nonlocal models is in the 
calculation of the angular diameter distances that determine 
E c and the relation between radial and angular scales. In par¬ 
ticular, the smaller values of H(z)/Hq in the Galileon and 
Nonlocal models lead to larger angular diameter distances, 
since these are oc / 1/H(z)dz. 

Figure[2]shows, for the same parameters as in Fig. [I] the ra¬ 
dial profiles of V 2 $ (solid) and (dashed) for the Galileon 
(upper left, red) and Nonlocal (upper right, blue) models. 
For both models, the figure also shows the predictions from 
models called QCDM (green), which are illustrative mod¬ 
els with the same background cosmology as their respective 
modified gravity models, but keeping the gravitational law of 
GR. Comparing the QCDM variants with the respective full 
models allows one to isolate the effects of the modified back- 
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ground from those due to the modified gravitational law. Since 
both QCDM models differ from ACDM only via the mod¬ 
ified background, their relative differences in the profiles of 
V 2< 1> ~ p and M(< R) are determined only by the 

background dependence of the density and mass profiles of 
the halos. As a result, the QCDM curves in Fig(2] follow the 
same behavior seen in the lower panel of Fig. [T] It is therefore 
more interesting to analyse the impact of the fifth forces. In 
the case of the Nonlocal model, we have seen that the modifi¬ 
cations to gravity can be parametrized by a scale-independent, 
but time-evolving effective gravitational strength, Gii n (t:), ac¬ 
cording to Eq. ( |34| ). Consequently, the effect of the fifth force 
in this model is to boost the amplitude of V 2< f> and $,a by the 
same amount on all scales. From Fig. |T] at z,i = 0.5 one has 
Gii n ~ 0.022, which corresponds to the difference between 
the predictions of the QCDM and Nonlocal models seen in 

Fig.0 

The effects of the fifth force in the Galileon model are 
slightly more complex due to the nonlinear screening mecha¬ 
nism. By comparing the predictions of the full Galileon model 
with those of the corresponding QCDM variant, one can iden¬ 
tify three regimes. The first is a "fully-screened" regime, R < 
O.IMpc/ft , where the effects of the fifth force are almost neg¬ 
ligible (<f> ss $ GR ), as seen by the overlap between the red and 
green sets of curves. The second regime is a "partly screened" 
regime which occurs on scales 0.1Mpc//i < R < 50Mpc/A 
On these scales, we can write $ = a(r)$GR^ where the func¬ 
tion a encapsulates the scale-dependence of the fifth force. 
Finally, on scales r > 50Mpc//i, the fifth force becomes com¬ 
pletely unscreened and one effectively has $ = Gii n (z)<f>GR, 
where Gy n is given by Eq. ( |23| > for the Galileon model. This 
translates into a constant boost in the values of V 2 <f> and '1>,« 
at large radii. The size of this boost in the unscreened regime 
of the Galileon model (15 — 20%) is larger than that in the 
Nonlocal gravity model (2—3%). This follows from the higher 
value of Gun in the Galileon model at z,i = 0.5, as shown in 
the middle planel of Fig. |T| 

The lower panels of Fig. [2] show the lensing convergence 
angular profiles for the Galileon (red) and Nonlocal (blue) 
models, as well as their respective QCDM (green) variants, 
plotted as the difference relative to ACDM and assuming 
z s = 1. These convergence angular profiles relate to the radial 
profiles of V 2 <1> by (i) the integration along the line of sight; 
(ii) the overall amplitude scaling set by £ c ; (iii) and also im¬ 
portantly, horizontal shifts caused by the fact that the same an¬ 
gular scales correspond to different distance scales at the clus¬ 
ter position because of the different D,j values. In the Galileon 
and Nonlocal model backgrounds, Dd becomes larger, and as 
a result, the same radial scales correspond to smaller angular 
scales. The net result of these effects is to reduce slightly the 
relative differences of k in the Galileon model w.r.t. ACDM, 
compared to the relative differences observed in V 2( I>. The 
same holds for the case of the Nonlocal model, for which the 
relative difference becomes also very weakly dependent on 
the angular scale (the slope of the curves is hardly noticeable 
in the scale of the figure). 



FIG. 3. Dependence of the factor T = Z M (zd, z s )/Z nd (zd, z s ), 
Eq. (|52j, on the source redshift z s for a ACDM model with f2 m o = 
0.4 (black), Galileon gravity (red) and Nonlocal gravity (blue). The 
cosmological background of these models is different from the fidu¬ 
cial ACDM model with D m o = 0.27 used by Ref. 11661 . In this 
figure, Zd, = 0.35, which is typical for the CLASH clusters. 

IV. METHODOLOGY 

We estimate cluster masses in Cubic Galileon and Nonlocal 
gravity cosmologies using the radially-binned lensing conver¬ 
gence profiles obtained from the reconstructions of the lensing 
potential for 19 X-ray selected galaxy clusters from CLASH 
ll65ll . In this section, we describe our methodology, paying 
particular attention to a number of subtleties that need to be 
accounted for to self-consistently compare the data with pre¬ 
dictions from the alternative models studied here. 


A. Cluster convergence profiles in alternative models 

We use the convergence profile data that was obtained for 
the CLASH clusters in Ref. Il66l . There, the analysis was 
performed with a numerical algorithm called SaWLens lf84l . 
which iteratively reconstructs the lensing potential for each 
cluster on a two-dimensional grid that covers the cluster field. 
The analysis is non-parametric, i.e., it makes no assumptions 
about the mass distribution of the cluster. We refer the reader 
to Refs. Il66ll84ll85l for the details about how SaWLens oper¬ 
ates. For the discussion here, what is important to note is that 
what SawLens actually reconstructs is the lensing potential 
scaled to a source redshift of infinity, t/>oo = ip(z s = 00 ), 
by assuming a fiducial cosmological model. We use to 
denote the lensing convergence associated with tpao, which is 
related to the convergence at the true source redshift, z s , via 

k Zs = Z m (z d ,z s )K oo, (47) 

where we use the subscript Zs to emphasize that n Zs corre¬ 
sponds to the convergence associated with z s . The function 
Z = Z(zd, z s ) transports the convergence from a source red¬ 
shift of infinity to the source redshift that corresponds to the 
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galaxies on each SaWLens grid cell/pixel (we use the words 
cell and pixel interchangeably). It is given by 


Z Rd (z d ,z s ) 


D^ooD 


fid 




(48) 


where the superscript fid indicates angular diameter distances 
that are calculated assuming the fiducial background cosmol¬ 
ogy and the subscript means that the calculation assumes 
that z s = oo. In the reconstruction process of Ref. 1661 . 
the fiducial cosmology is a ACDM model with Q m o = 0.27. 
From hereon, we use to denote the convergence profiles 
obtained in this way, where the superscript fid makes it explicit 
that the data is linked to the fiducial model. It is therefore im¬ 
portant to investigate the extent to which the profiles can 
be used in studies of alternative cosmologies. 

Consider the case that we wish to estimate the lensing 
masses of the CLASH clusters in a model with a cosmological 
background that is different from the fiducial model originally 
used to analyse the observations in Ref. l66l . In principle, we 
could suitably modify the SaWLens algorithm to reconstruct 
the convergence maps in the alternative model, K'lf, instead of 
However, this would not be practical as it would imply 
rerunning the entire analysis pipeline for different background 
cosmologies. A more economical strategy is to note that the 
two convergence maps, rff and can be related by 


Z dd (z d ,z s )^ = Z^(z d ,z s ) K : 


alt 
v oo 1 


(49) 


where Z alt is defined as in Eq. ( |48| ) but with the distances 
calculated in any alternative, and not the fiducial, cosmology. 
The above equation holds (up to a correction that we discuss 
in the next subsection) since both Z dd n^ and Z alt K^ corre- 

f ond to K Za , i.e., the convergence at the true source redshift 
Using the above equation, the radially binned convergence 
profiles obtained using the fiducial cosmology in Ref. (66), 
K oo W> can b e directly compared to the prediction of the al¬ 
ternative model k^‘, provided the latter is multiplied by the 
factor Z aXt {z d , z s )/Z Rd (z d , z s ). This is the approach that we 
adopt in this paper. Specifically, we aim to obtain constraints 
on M 2 oo, C 200 and 9 ra0 in non-fiducial backgrounds by mini¬ 
mizing the x 2 quantity 


X 2 = (50) 


where 


Vi — — T^(M 2O o,C2oo,fi m o,0i), (51) 

is the /-th entry of the vector V ; i is the reconstructed 
lensing convergence in the i -th radial bin, 0 ,; C K is the covari- 


2 For example, for ACDM with an alternat ive background and for fixed 
surface mass density, for simplicity, Eq. |49| becomes E(0)/Ejl 14 = 


E(0)/E 


fid 


■ EJ“ = E; 


fid 


ance matrix of the radially binned data 0 and for brevity of 
notation, we introduce the scaling factor 


T C(z d ,z s ) 


Z alt (z d ,z s ) 
Z dd (z d ,z s )' 


(52) 


In Eq. (5l) , is given by Eq. (l2) for ACDM, Eq. (25 1 for 
Galileon gravity and Eq. (36) for Nonlocal gravity, but using 
z s = 00 , in the calculation of £ c . 

Unless otherwise specified, we assume flat priors on the 
free parameters, U m0 € [0.1,0.5], M 20 0 € [0.3,3.0] x 
10 15 M 0 /h and c 20 o € [1,8]. 


B. The validity of Eq. (49) and the choice of source redshifts 


As discussed above, is reconstructed by applying the 
transformation of Eq. (47) in each cell of the SaWLens grid 
that covers the cluster field. In this process, the value of z s 
is determined by the redshift of the galaxies used to measure 
the ellipticity field at that pixel, or by the redshift of the galax¬ 
ies associated with a given multiple image system. On the 
other hand, our methodology is based on Eq. (49) , in which 
one scales the lensing quantities from z s = 00 to a source 
redshift z s , but neglects the redshift distribution of the back¬ 


ground lensed galaxies. The validity of Eq. (49 1 then becomes 
linked to the impact of the spread of the redshift distribution of 
the source galaxies across each cluster field. For the CLASH 
clusters analysed in Ref. (66ll . the redshift distribution of the 
background galaxies is manifest in four main aspects: 

(i) in the weak lensing regime, different ellipticity pixels 
are associated with different source redshifts since the shapes 
are measured using different galaxies across the cluster field; 

(ii) related to the above, the ellipticity of each pixel results 
from a local average of neighbouring galaxy shapes, which 
can have different redshifts; 

(iii) the ellipticity field used by SaWLens is a combined 
catalog of measurements from space- and ground-based tele¬ 
scopes, which probe different galaxy redshift ranges. The 
measurements of these two catalogs (see Ref. (66l ) are cor¬ 
rected for this, but assuming the fiducial cosmology; 

(iv) in the strong lensing regime, each pixel is associated 
with the redshift of the multiple images contained within it, 
which can be different in different multiple image systems for 
the same cluster and also different from the galaxy popula¬ 
tions used in the weak lensing measurements. 

To get a feeling for the size of our approximation, we show 
in Fig. [3] the z s dependence of the factor T (Eq. (52) ) for 
the Galileon (red) and Nonlocal (blue) models, and a ACDM 
model with = 0-4 (black). The quantity T encapsulates 


3 The bootstrap realizations used to derive the covariance matrices in 
Ref. (66) also make use of the fiducial cosmological background. Here, 
we use the errors as obtained for the fiducial cosmology and do not attempt 
to estimate the dependence of the covariance matrix on the assumed cos¬ 
mology. This does not alter our conclusions as this choice only affects the 
precise size of the confidence intervals, without introducing any important 
systematics. 
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all of the dependence on z s in the x' 2 minimization used to 
estimate the cluster parameters. For illustrative purposes, we 
choose Zd = 0.35. This corresponds roughly to the mean red- 
shift of the CLASH clusters (cf. Table [I]), although the exact 
value is not important for the discussion here. We note that 
what is relevant is the slope of the curves and not their abso¬ 
lute value. Consider for the sake of argument an extreme case 
where the source galaxies are distributed between z s = [1,3], 
but that we choose to use z s = 2 in Eq. ( [52] >. Focusing on the 
case of the Galileon model, we have that T(z s = 1) sa 1.019, 
T(z s = 2) « 1.011 and T (z a = 3) « 1.007. These val¬ 
ues differ by no more than w 1%, and hence our choice of 
z s should not lead to serious biases in the results. The error 
would be even smaller in the Nonlocal model or ACDM with 
n m0 = 0.40, since in these cases the Y(z s ) curves are shal¬ 
lower than in the Galileon case. The error of neglecting the 
redshift distribution becomes smaller for higher values of z s , 
for which the curves in Fig. [3] become visibly flatter. This is 
relevant for strongly lensed systems, which tend to be associ¬ 
ated with galaxies at higher redshifts. 

In cluster weak lensing studies, it is common to determine 
an effective source galaxy redshift, defined as 


Dd s . , _ / Dds \ 


(53) 


where ( Dd s /D s ) is an average over all source galaxies. Ref¬ 
erence (66l quotes z s , e s values for the CLASH clusters (see 
also Ref. (67)). For example, Abell 209 (z d = 0.206) has 
( Dds/D s ) = 0.75 ± 0.04 (lu), which corresponds to z Sje ff = 
1.03^o'i5 (this estimate comes from Table 3 of Ref. (671). 
This uncertainty on z Sie ff is much smaller than our rather ex¬ 
treme example above (z s = 2 ± 1), which further convinces 


us that the approximation of Eq. (49 1 is a good one. For com¬ 
pleteness, we note that the determination of these values of 
z s ,efi involves knowledge of the background cosmology, and 
hence they are also model dependent. However, again tak¬ 
ing Abell 209 as an example, in the Galileon model one has 
Dds/D s (z S:e ff = 1.03) = 0.76, which is well within the 
uncertainty (±0.04) quoted above for this cluster. We can 
therefore neglect this model dependency and use the values of 
z v ,ff listed in Ref. (66). In particular, in our \ 2 minimization, 
we shall use the effective source redshift values found for the 
background sources of the ground-based ellipticity measure¬ 
ments, which we list in Table [I] 

To summarize this discussion, although Eq. ( [49] ) is only ap¬ 
proximate, the results shown in Fig. [3] suggest that our results 
are insensitive to the exact choice of z s . 


C. Other subtleties in using cluster lensing data to test gravity 

Before proceeding further into estimating the CLASH clus¬ 
ter masses in modified gravity models, we discuss some other 
subtle issues that may arise when combining current lensing 
modelling techniques with modified gravity. Although it turns 
out these other issues do not play a direct role in the results 
of this paper, we believe such a discussion is instructive and 


leads to a clearer and broader understanding of the results of 
this and other work in the literature. 


1. Parametric vs. nonparametric analysis 

The non-parametric reconstruction of the lensing potential 
used in this paper builds solely upon the observed lensing con¬ 
straints, without making any assumptions about the mass dis¬ 
tribution of the cluster. Such a model-independent [^] method 
is particularly well suited to modified gravity studies. Con¬ 
sider the alternative scenario of a parametric approach. In this 
case one starts by making an Ansatz about the mass distribu¬ 
tion in the cluster. Typically, this can involve describing the 
main dark matter distribution using a single (or more in the 
case of mergers) NFW profile. Then, one could also model 
substructure by identifying the position of the most massive 
cluster galaxies and assigning them a given density profile, 
(see e.g. Refs. lf86Ul01lD . The free parameters of such a mass 
model are then iterated over until the lensing constraints are 
satisfied. In the context of modified gravity there are at least 
two subtle issues associated with such a parametric lensing 
analysis. First, in order to compute the lensing effects due 
to the postulated mass distribution one must assume a theory 
of gravity: for the same mass distribution, different models 
of gravity could induce different lensing effects. Parametric 
methods are therefore biased towards the assumed theory of 
gravity. Second, the lensing properties of a given point in the 
cluster field are determined by the sum of the lensing signal 
predicted by each element of the mass model (main halo plus 
the substructures). This superposition is valid in GR (which 
is linear in the Newtonian limit), but not necessarily in alter¬ 
native (typically nonlinear) models of gravity. These issues 
can be circumvented if one reconstructs directly the lensing 
potential and its derivatives but not the mass distribution. It is 
for this reason that we choose to use the SaWLens results of 
Ref. (66) in our analysis. 

For clarity, it is worth pointing out that the problem with 
parametric mass modelling is less severe if it is only applied 
to the strong lensing part of the data analysis. In this regime 
(well within /Goo), the effects of the modifications to grav¬ 
ity in a model like the Galileon may be small by virtue of 
the screening mechanism. As a result, the assumptions made 
in the data analysis may turn out to be a good approxima¬ 
tion. This is the case for some of the Galileon-type models 
explored in Refs. 1 102 , 11031 . which were compared to (para¬ 
metric) cluster data from Refs. [f97] , |98l . However, it is rea¬ 
sonable to expect that the screening efficiency is different in 
different models of gravity. For instance, the Nonlocal grav¬ 
ity model predicts a constant enhancement of the gravitational 
strength on all relevant scales. One must therefore be always 
cautious and check for the compatibility of the data analysis 
with the specific theory of gravity to be tested. 


4 Apart from the issue of the fiducial cosmological background model dis¬ 
cussed above. 
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2. Interpretation of stacked cluster leasing profiles 

To overcome systematic effects due to intervening struc¬ 
ture, cluster substructure and cluster asphericity, it has be¬ 
come common to build average (stacked) lensing profiles by 
using cluster lensing data from independent lines of sight 
l97l 198 . 1041 . The averaged profiles are then fitted again to 
infer an average mass and concentration that characterizes the 
stack. From a conceptual point of view, the same procedure 
can be applied assuming modified gravity models. Here, we 
comment that the interpretation of the stacked data may be 
somewhat more complex due to the effects of modified grav¬ 
ity. Consider for simplicity the stacking of the convergence 
radial profiles of N clusters at redshifts Z\ jy with mass and 
concentration values Mi. jv and ci. jv, respectively. The back¬ 
ground galaxies can be assumed to lie at the same source red- 
shift. For instance. Ref. ED stacks four massive clusters by 
co-adding (with some weighting) their profiles. The resulting 
mean profile is then refitted to determine a mean mass and 
concentration of the stack. Now consider fitting such a stack 
to two gravity models which display different time evolution 
for an unscreened gravitational strength. For these two mod¬ 
els, clusters located at different redshifts would contribute dif¬ 
ferently to the mean mass/concentration estimate since their 
lensing signal is amplified differently. For such a scenario, 
an interesting analysis would be to split the stack into smaller 
ones binned by cluster redshift and check for differences in 
the resulting mean mass/concentration of the smaller stacks. 
The situation becomes even more complex (but interesting) 
in models with screening, due to its scale-dependence, whose 
efficiency is in general redshift dependent as well. 

We stress that the above issues do not pose a serious prob¬ 
lem to using stacked data to test modified gravity, but simply 
that the extra physics can enrich the interpretation of the re¬ 
sults. In this paper, however, we shall not be concerned with 
these issues since we fit each of the CLASH clusters individ¬ 
ually. 


V. RESULTS: LENSING MASS ESTIMATES 

In this section, we present our main results for the mass and 
concentration estimates of the CLASH clusters in the three 
models of gravity we consider. We also discuss the impact of 
our results in the context of recently proposed observational 
methods to test gravity on large scales. 


A. The impact of H m0 on cluster lensing mass estimates 

Figure [4] shows the constraints in the C 200 -M 200 (left), 
fl m o-M 2 oo (middle) and fl m o-c 2 oo (right) planes for the RX 
J2129 (blue) and CL J1226 (red) CLASH clusters, assuming 
a ACDM model. The solid contours show the parameter con¬ 
straints obtained using a three-dimensional grid search. The 
gray dots correspond to the points accepted in the chains of 
a Monte Carlo exploration of the parameter space using the 


recently developed and publicly-available Cosmo SI ^pack¬ 
age 11051 . CosmoSIS is a highly modular parameter estima¬ 
tion code, to which we have added likelihood modules based 
on the convergence profiles of each of the 19 CLASH clus¬ 
ters. Given that our grid code and the MCMC sampler in 
CosmoSIS are independent ways of sampling the parame¬ 
ter space, it is reassuring that their results agree as shown in 
Fig. 0 We will be making the CLASH likelihood modules 
developed in this paper publicly available in future releases of 
the CosmoSIS standard library. 

The middle and right panels of Fig. [4] show that the data 
from each cluster does not place any meaningful constraints 
on fi m o (this conclusion holds for all of the other CLASH 
clusters as well). Furthermore, there is no clear degeneracy 
between fi m0 and the cluster parameters M 2 00 , c 2 00 , which 
indicates that the constraints on the mass and concentration 
are unbiased w.r.t. a particular choice of fi m 0 . This point is 
also illustrated by the dashed contours in the left panel of 
Fig. 0 which show the constraints when the cosmological 
matter density is held fixed at fl m 0 = 0.27. The compari¬ 
son of the solid and dashed contours shows that there is no 
significant deterioration of the constraints when one allows 
fl m 0 to be a free parameter. The cases of CL J1226 (red) 
and MACS J0744 (not shown) are those for which the dete¬ 
rioration is the most pronounced (although still very small). 
For the rest of the clusters, the two sets of contours (solid 
and dash) are barely distinguishable, just as in the case of RX 
J2129 (blue). This illustrates that, in observational determi¬ 
nations of cluster masses from lensing, the impact of assuming 
a specific value of fl m0 (as is often done in the literature) is 
negligible [^] 

To help understand the above result (and others that will 
follow), it is instructive to look at the impact of M 2 00 , c 2 oo 
and fl m0 on the lensing convergence profiles. This is shown in 
Fig.[5]for ACDM (black), Galileon gravity (red) and Nonlocal 
gravity (blue). For all models: 

• an increase in c 2 oo (keeping all other parameters fixed) tilts 
the convergence profile, boosting its amplitude in the inner 
regions and lowering it in the outer regions of the halo; 

• larger mass values enhance the lensing convergence on 
all scales shown, with the effect being slightly more pro¬ 
nounced at large radi0 


5 The bitbucket webpage of the CosmoSIS code is 
https://bitbucket.org/joezuntz/cosmosis/wiki/Home . 

6 The changes in the constraints when one allows o to vary are so small 
that they are of the same magnitude as the error of the approximation of 
Eq. J49k which is expected to be negligible in any case, as discussed in 
Sec.pVBl 

7 Note that for 0 > 10 2 arcsec, the effects of increasing halo mass start to 
become degenerate with the effects of decreasing the concentration. Strong 
lensing analysis, which probes the inner most regions of the cluster, there¬ 
fore helps to break this concentration-mass degeneracy. Reference ED 
also determines the lensing profiles of some of the CLASH clusters used 
in Ref. 1661 . However, the latter includes strong-lensing constraints in the 
analysis, which allows for more accurate concentration estimations. This 
is why we choose to work with the results of Ref. 166 ). 
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FIG. 4. Two-dimensional marginalized constraints on the C200 — M 200 (left), H m o — M 200 (middle) and H m o — C200 (right) planes, for the 
CLASH clusters RX J2129 (blue) and CL J1226 (red), assuming a ACDM cosmology. The solid contours depict the marginalized 68% and 
95% confidence limits when all three parameters are varied. The dashed contours in the left panel are the same as the solid ones, but with the 
cosmological matter density held fixed at fl m 0 = 0.27 (barely noticeable in the case of RX J2129). The gray dots show the points accepted in 
the Monte Carlo chains built by the CosmoSIS code. 
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FIG. 5. Lensing convergence k(9) for ACDM (black), Galileon 
gravity (red) and Nonlocal gravity (blue) models, plotted as the 
difference relative to their respective models with base parameters 
M 200 = 5 x 10 14 M 0 //t, C200 = 4 and fi m o = 0.27. The blue 
and black curves are almost overlapping. We assume Zd = 0.5 and 
z s = 1.0. The dashed, dotted and dot-dashed lines correspond, re¬ 
spectively, to cases for which the values of A/200, C200 and are 
enhanced by 10%. To guide the eye, the shaded region indicates the 
radial scales beyond R 200 in the base ACDM model. 


• A boost in 1 l rrM has three main effects. First, it increases 
the value of p c (z ) at the cluster redshift, which effectively 
increases the density of the halo for fixed M 200 (recall the 
discussion of Sec. m- This should boost the convergence at 
small radii. Second, an increase in fl m o also has an impact 
on the cosmological distances that enter £ c , which has the 
overall effect of lowering the amplitude of k, by a constant 
factor. Third, a change in f2 m0 also changes the angular di¬ 
ameter distance to the cluster, Dd, which causes horizontal 


shifts in the lensing convergence, when plotted as a function 
of the observed angular scales 6 = Dd/r. 

The net effect of varying f2 m o results in almost no visible 
shift (dot-dashed) in the amplitude of k for the radial scales 
that are better probed by the CLASH radially binned profiles, 
9 < 500 — 700 arcsec (cf. Fig. [7}. This is why the constraints 
on fl m0 in Fig. 0 are so weak. We note also that since the ef¬ 
fect of fl m g is always subdominant w.r.t. the effects of varying 
M 200 and C 200 (in terms of fractional change), it is unlikely 
that future lensing data from experiments such as EUCLID 
Iil06l and LSST 111071 will change this conclusion. 

From hereonin, we fix the cosmological matter density to be 
Q m0 = 0.27 for all our models. In this way, the ACDM model 
becomes the fiducial one used in Ref. l66l . This value is also 
consistent with the CMB observational constraints for the Cu¬ 
bic Galileon and Nonlocal models as found, respectively, by 
Ref. 10 and Ref. (28). 


B. Cluster lensing masses in Galileon and Nonlocal gravity 

Figure [6] shows the constraints on the r^oo-AGoo plane ob¬ 
tained for each of the CLASH clusters in ACDM (black), 
Galileon (red) and Nonlocal gravity (blue) cosmologies. The 
dots indicate the position of the best-fitting values. The best- 
fitting lensing convergence profiles are shown in Fig.[7](what 
is shown is Tk oo ( 0)). The concentration-mass relation of the 
CLASH clusters for the three models is shown in Fig. |8] to¬ 
gether with results from N-body simulations ll29l |61] Tl 081 . 
First, we note that our cluster mass and concentration esti¬ 
mates for ACDM are in agreement with those obtained in 
Ref. l66l . Second, these three figures all show that the con¬ 
straints on the cluster parameters are, within errorbars, the 
same in the three cosmological models. Although there are 
tiny differences in the resulting best-fitting values of M 200 
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FIG. 6. Two-dimensional 68% and 95% confidence limits on the C 200 — M 200 plane for all of the CLASH clusters assuming ACDM (black), 
Galileon gravity (red) and Nonlocal gravity (blue). The position of the best-fitting points is marked by the dots, and their respective y 2 values 
are shown in each panel. 


and C 200 for the three models (< 5%), they all lie well within 
the 1 a limits (whose precision varies within ~ 50% — 80%). 
The shapes of the contours are also remarkably similar and the 
goodness-of-fit is essentially the same in all models, as can be 
seen by comparing the respective \ 2 values in Fig. [6] In Fig. [7] 
one notes that for almost all of the clusters, the best-fitting 
convergence profiles underpredict the data points at large an¬ 
gular scales (although well within the errorbars). However, 
close to the edge of the clusters, the contribution from the 
surrounding large scale structure may have a non-negligible 
impact. This can partly explain why the data points tend to go 
up at large scales, as investigated, for instance, in Ref. 111091 . 

The shaded bands in Fig. [8] show the best-fitting mean 
concentration-mass relations found in N-body simulations for 
the ACDM (gray) model in Ref. 111081 the Cubic Galileon 
model (red) in Ref. 11611 and the Nonlocal model (blue) in 
Ref. |29l . In these bands, the lower and upper bounds cor¬ 
respond, respectively, to the relations at z = 0.666 (a = 0.60) 
and z = 0 (a = 1) (this redshift range is approximately that 
of the CLASH clusters). Figure [8] shows that there is good 
agreement between the simulation results and the concentra¬ 
tion/mass estimates of the CLASH clusters in the three models 
of gravity. However, there are a number of issues that prevent 
a direct comparison between the simulation results and the 


See Fig. 9 of Ref. [66} for the comparison of the CLASH C 200 -M 200 rela¬ 
tion in ACDM with other relations in the literature. 


estimated concentration and mass values. First, the shaded 
bands of the Galileon and Nonlocal models in Fig.[8]have been 
extrapolated to masses larger than the mass range used to fit 
the best-fitting concentration-mass relations in the simulations 
of Refs. ED and ll29l . Second, the concentration-mass rela¬ 
tion was fitted using all haloes, without applying any selec¬ 
tion criteria to consider only relaxed ones mm This may 
be particularly relevant for the CLASH clusters, which are 
characterized by regular X-ray surface brightness morpholo¬ 
gies l65l . and are therefore expected to be relaxed and close 
to virial equilibrium (see also Refs. Ifim HT2l for a recent 
discussion on the impact of baryonic processes in the den¬ 
sity profiles of clusters). Third, the concentration-mass rela¬ 
tion in the simulations was obtained by fitting NFW profiles to 
the three-dimensional spherically averaged mass distribution 
of the haloes, whereas the symbols in Fig. [8] are the values 
obtained by also assuming spherical symmetry, but fitting to 
two-dimensional (projected) lensing convergence profiles (see 
e.g. Sec.6.2 of Ref. l66l for an analysis of the impact of this 
projection bias in the CLASH sample). Finally, the upper and 
lower bounds of the bands correspond to the mean relation 
found in the simulations, but the intrinsic scatter around the 
mean concentration-mass relation should also be taken into 
account. Nevertheless, to guide the eye, we opted to keep the 
simulation results in Fig. [8] but advise that further work is 
needed before performing a more thorough comparison (see 
the analysis of Ref. 1113l in ACDM models for an illustration 
of the steps to follow). 
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FIG. 7. Best-fitting lensing convergence profiles, k(6) = Tk.oo, obtained for all of the CLASH clusters assuming ACDM (black), Galileon 
gravity (red) and Nonlocal gravity (blue). The green dots are the radially binned data as described in Ref. [66ij and the errorbars are the square 
root of the diagonal entries of the covariance matrix of the data. To guide the eye, the dotted vertical lines indicate the inferred values of R 200 , 
which are barely distinguishable for the three models. 


The left panel of Fig. [9] shows the best-fitting lensing con¬ 
vergence for all of the CLASH clusters in the Galileon (red) 
and Nonlocal (blue) cosmologies, plotted as the respective dif¬ 
ference to the best-fitting profiles in ACDM. As expected 
from the above results, on the scales that are probed by the 
CLASH data, 6 < 500 — 700 arcsec, the three models are 
in very good agreement. In the case of the Galileon model, 
this is because the screening is very effective on these scales 
inside f?2oo (Ref. B1 141 finds a similar screening efficiency in¬ 
side i ?200 f° r DGP gravity, which employs also the Vainshtein 
mechanism). This can be noted by comparing the enhance¬ 
ment in the amplitude of k on larger scales, where the screen¬ 
ing becomes less efficient. In the case of the Nonlocal model, 
although the modifications to the gravitational strength are not 
screened, they are not strong enough to have a significant im¬ 
pact on the lensing convergence profiles. We therefore con¬ 
clude that, for the case of the CLASH clusters analysed here, 
the impact of modifying the lensing gravitational potential ac¬ 
cording to Cubic Galileon or Nonlocal gravity is completely 
negligible in the estimation of their lensing masses. 

Before this paper, there have been other works investigat¬ 
ing the impact of Galileon-like effects on cluster lensing pro¬ 
files. Reference B1Q3I used a parametrization of the fifth 
force, which was constrained using the stacked cluster lens¬ 
ing shear profiles from Refs. £2112D . Their parametrization 
encompasses the Galileon model studied here, and for this 
case, there is good agreement with our conclusions. Refer¬ 


ence B1021I performed similar investigations, but took as a test 
case a model inspired by massive gravity. The authors do not 
consider the time evolution of the cosmological background 
and their equations of motion include higher order terms than 
the Galileon model studied here, which prevents a direct com¬ 
parison with our results. A particularly interesting feature de¬ 
scribed in Ref. 11021 is that some model parameters predict a 
"dip" in the amplitude of the convergence profiles, which hap¬ 
pens for r < .R200 ( see e.g. Fig. 4 of Ref. 111021 ). These scales 
are sufficiently well probed by the lensing data, which allows 
some of these specific models to be ruled out already. Such 
features, however, do not show up in our convergence profiles 
for the Cubic Galileon model (cf. Fig.[7]l, but which is a differ¬ 
ent model from the one studied in Ref jfl02ll . In the context of 
f(R) models, although the lensing signal is not modified di¬ 
rectly, it can be via modified mass distributions. For instance, 
the enhanced dynamical potential in these models boosts the 
accretion rate of matter onto the clusters, which results in an 
excess of mass in their infall region. This can be probed with 
lensing measurements, as was done in Ref. mi. Since the 
Galileon and Nonlocal models also directly modify the dy¬ 
namical potential, then in principle, similar investigations can 
be performed. In the case of these models, the lensing signal 
should be amplified both by the excess of mass that surrounds 
the cluster, and by the intrinsically enhanced lensing effects. 
A detailed investigation of this is, however, beyond the scope 
of the present paper. 
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FIG. 8. Concentration-mass relation of the CLASH clusters assum¬ 
ing ACDM (black), Galileon gravity (red) and Nonlocal gravity 
(blue). The errorbars indicate the marginalized 68% confidence lim¬ 
its. We use different symbols for different clusters to facilitate the 
identification of which cluster is which across the three models. The 
shaded bands indicate the mean concentration-mass relations from 
N-body simulations between 2 = 0.66 (lower bound) and 2 = 0 
(upper bound) found for ACDM (gray) in Ref. lll08l . Cubic Galileon 
model (red) in Ref. ED and Nonlocal model (blue) in Ref. CH. 


C. The connection with tests of gravity 


1. Dynamical masses from the phase-space density around 
massive clusters 


Recently, the authors of Refs. 11611191 have proposed meth¬ 
ods to test the law of gravity on Mpc scales by using informa¬ 
tion from the galaxy velocity field in the infall regions around 
massive clusters (see also Ref. 11161 ). These techniques were 
designed with models of gravity that modify the dynamical 
potential (i.e. that felt by nonrelativistic objects like galaxies), 
but do not modify the lensing potential (i.e. that felt by rela¬ 
tivistic particles like photons). Popular models such as f(R) 
and DGP gravity fall in the above category, and as such, the 
lensing mass estimates, M\ en , for these models would auto¬ 
matically be the same as in GR. On the other hand, the ve¬ 
locity dispersion of surrounding galaxies as they fall towards 
the clusters would be affected by the modifications to gravity. 
Therefore, if one would interpret these observations assuming 
GR, then one would infer dynamical masses, Md yn , which 
are different from those estimated using lensing. A mismatch 
in the estimates of the lensing and dynamical masses would 
therefore be a smoking gun for modified gravity lfT3l IT5 .. fTTI 
(see, however. Ref. Ill 181 for a discussion of how complica¬ 
tions associated with assembly bias could affect these tests). 

The merit of the test of gravity described above becomes 
less clear when applied to models that also modify the lens¬ 
ing potential. Consider, for simplicity, a model that boosts the 
dynamical and lensing potential by the same constant factor, 
a > 1, i.e. ‘f’dyn = ®icn ~ a$ GR . In such a model, the 
mass of a cluster inferred from the surrounding galaxy veloc¬ 


ity field would be biased low w.r.t. GR. This is because, due 
to the enhanced gravitational strength felt by the galaxies, the 
cluster does not need to be as massive as in GR to accelerate 
the galaxies by the same amount. Following the same reason¬ 
ing, the lensing mass estimates would also tend to be biased 
low compared to GR: due to the fifth force felt by the photons, 
the cluster can be less massive to induce lensing effects of the 
same magnitude. In such a model, both Md yn and M\ en shift 
in the same direction. This therefore makes it harder to tell 
the two values apart and hence, harder to detect a signature of 
modified gravity. 

The Galileon model also modifies the lensing potential, but 
adds complexity to the case described above in the sense that 
the modifications to gravity are scale dependent with screen¬ 
ing inside the cluster radius. Just outside the cluster radius, 
the screening becomes less efficient and the fifth force sig¬ 
nificantly boosts the lensing convergence, as shown in Fig. 
9. Although these larger scales are not accurately probed by 
the current cluster lensing data, they correspond roughly to 
the regions associated with galaxy infall, 2Mpc jh < r < 
20Mpc /h. For these radial scales, the right panel of Fig. 9 
shows that the total force profile which surrounds the CLASH 
clusters in a Galileon cosmology can be up to 10% — 40% 
higher than in ACDM. As a result, galaxies located at these 
distances from the cluster center should feel the boost in the 
total force, which should translate into their velocity distri¬ 
bution. On these scales, both the lensing and dynamical 
masses would be different in a Galileon cosmology compared 
to ACDM. Inside the cluster radius, on the scales that are 
probed by the CLASH data, 6 < 500 — 700 arcsec, the left 
panel of Fig. 9 shows that the differences in the convergence 
profiles compared to those in ACDM are small enough for 
the mass estimates to be almost the same in the two mod¬ 
els. Therefore, inside the cluster radius, this leaves us with 
a similar picture to that in f(R) or DGP models: the lensing 
mass estimates are not affected by the modifications to grav¬ 
ity, but dynamical mass estimates using infalling galaxies are 
changed. We therefore conclude that, despite it being a model 
that modifies the lensing potential, the fact that dynamical and 
lensing mass estimates are sensitive to radial scales of differ¬ 
ent screening efficiency allows the Cubic Galileon model to be 
tested by the methods proposed in Refs. EMU. 

In the case of the Nonlocal model, although the lensing 
mass estimates are also practically the same as in ACDM, the 
enhancement of the force profile on scales 2Mpc /h < r < 
20Mpc/ft is kept below the ~ 5% level. This makes it more 
challenging for this model to be tested by these methods. 


2. Galaxy-galaxy lensing 

The left panel of Fig. [9] also shows that although the con¬ 
vergence profiles are very close in the three models for R < 
f? 200 > they can be visibly higher (by ~ 20 — 80%) in the 
Galileon model on larger scales. The enhanced lensing sig¬ 
nal outside dark matter haloes in Galileon-like models has 
been analysed by Refs. El SSI, but in the context of theo¬ 
ries that emerge from massive gravity scenarios |[33ll38l . In 
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FIG. 9. (Left) Best-fitting lensing convergence profiles, k(9) = Tkoo, for the all of the CLASH clusters in the Galileon (red) and Nonlocal 
(blue) gravity models, plotted as the difference relative to the best-fitting profiles in ACDM (black). To guide the eye, the shaded band 
represents approximately the regions that lie beyond R 200 for all clusters. (Right) Same as the left panel, but for the total force profile 
The shaded band encloses the scales 2Mpc /h < R < 20Mpc//i which are approximately those associated with the infall of surrounding 
galaxies. 


particular, the authors investigate the possibility of such a 
signal being detected in galaxy-galaxy lensing observations 
(see e.g. Ref. Il47l ). The latter can be measured by cross- 
correlating the position of foreground galaxies (the lenses) 
with their background shear field. Our results in Fig. [9] are 
in good qualitative agreement with the solutions explored in 
Refs. 14511461. For instance, we also find the appearence of a 
bump in the relative difference to ACDM, which we checked 
occurs at ~ 10/Con- Quantitatively, the comparisons become 
less straightforward. On the one hand, in this paper we show 
the results for cluster mass scales between ss [0.5,1.5] x 
10 15 Mq//i, which are higher than the galaxy group mass 
scales (10 13 — 10 14 M 0 //i) probed in Refs. J45![4(3- More¬ 
over, our models also differ at the level of the cosmological 
background, exact screening efficiency and time evolution of 
the linearized effective gravitational strength. Nevertheless, it 
seems reasonable to expect that the predictions of the Galileon 
model studied here are also likely to be scrutinized by galaxy- 
galaxy lensing observations. A more detailed investigation of 
the model predictions for galaxy-galaxy lensing is beyond the 
scope of the present paper. 

In the case of Nonlocal gravity, the modifications to the 
lensing convergence are small (< 5%) on all scales, which 
makes it much harder to distinguish from standard ACDM 
with galaxy-galaxy lensing data. 


3. Weak lensing on larger scales 

The picture depicted in the left panel of Fig|9]that the lens¬ 
ing signal gets significantly enhanced on larger scales in the 
Galileon model should, in principle, also have an impact on 
the lensing of CMB photons. Indeed, Refs. BED have shown 


that the amplitude of the CMB lensing potential angular power 
spectrum, Cf^, is very sensitive to the modifications to grav¬ 
ity in the Galileon model. To the best of our knowledge, the 
effect of the Nonlocal model on the CMB lensing potential 
power spectrum has never been investigated in detail. How¬ 
ever, since the modifications to gravity on large scales are 
not as strong as in the Galileon model (cf. middle panel of 
Fig. [lj), the effects on the amplitude of Cf^ should be less 
pronounced. 

By the same reasoning the weak lensing cosmic shear 
power spectrum should also be sensitive to the modifications 
to gravity in the Galileon model, but less so in the Nonlocal 
case. Again to the best of our knowledge, cosmic shear data, 
such as that gathered by the CFHTLens Survey 111 191 , have 
never been used in direct tests of the models studied here, al¬ 
though Refs. [48, 49, 1201 have used these data to constrain 
general parametrizations of modified gravity. 


VI. SUMMARY & OUTLOOK 

We have estimated cluster lensing masses in alternative the¬ 
ories of gravity that modify the lensing gravitational potential. 
For this, we varied the mass (M 2 00 ) and concentration (C 200 ) 
of NFW haloes to fit the predicted lensing signal in modified 
gravity to the radially binned lensing convergence profiles ob¬ 
tained from non-parametric reconstructions of the lensing po¬ 
tential for 19 X-ray selected clusters from the CLASH survey 
EMO. The methodology we adopted is similar to that first 
employed in Ref. Il66l in the context of GR. 

We focused on the Cubic Galileon and Nonlocal models, 
which modify the gravitational law in qualitatively different 
ways. In the case of the Nonlocal model, the modifications 
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to gravity can be parametrized by an effective time varying 
gravitational strength, which is independent of length scale 
(cf. Eq. (|35]) and Fig. [TJ. In the Galileon model, the gravi¬ 
tational law can also be parametrized by a scale-independent 
gravitational strength on large scales (cf. Eq. ( |23| and Fig. [TJ. 
However, close to massive bodies, the Vainshtein mechanism 
(manifest in nonlinearities in the equations of the model) 
introduces a scale dependency to the gravitational strength, 
which acts to suppress the amplitude of the modifications 
w.r.t. standard GR/ACDM. The cosmological background in 
both models is also modified relative to ACDM (cf. Fig. [TJ, 
which has an impact on the conversion between angular and 
radial scales, p c {z) and also on the values of £ c . 

We paid particular attention to the compatibility of the data 
analysis with the modified gravity models we wished to test. 
Namely, we pointed out that the CFASH cluster convergence 
profiles obtained by Ref. Il66l are particularly suited for mod¬ 
ified gravity studies since the analysis makes no a priori as¬ 
sumptions about the mass distribution, and it is therefore less 
model dependent. If one constructs first a mass model for the 
cluster, then one must postulate a theory of gravity to compute 
the lensing signal associated with that mass distribution. The 
lensing convergence maps obtained from such an approach 
could therefore be biased towards the assumed theory of grav¬ 
ity, which would prevent a direct comparison with other mod¬ 
els. We have also pointed out that the analysis of Ref. l66l is, 
however, not completely model independent, as it assumes a 
fiducial ACDM background model (fl m o = 0.27) to compute 
angular diameter distances. In Secs. |IV A| we explained that 
this extra model dependency can nevertheless be taken into 
account by applying a correction factor, Eq. (52 1 , to the con¬ 
vergence profiles predicted by models with different cosmo¬ 
logical backgrounds. This correction factor holds under the 
approximation that all background source galaxies lie at the 
same redshift, which as we argued in Sec. IV B turns out to 


be a good approximation with negligible impact on our con¬ 
clusions. 

Our main results can be summarized as follows: 


• Although (!,,,() is a parameter that enters the calculation 
of the lensing convergence, its impact is very small com¬ 
pared to the size of the effects of M 200 an d C 200 (cf. Fig. [5j. 
This means that assuming a particular value for fl m 0 does 
not introduce any significant biases in the cluster mass and 
concentration estimates. We have shown this explicitly for 
ACDM by simultaneously varying H, m o, M 200 and C 2 oo> 
and found barely any difference from the constraints on the 
cluster parameters obtained when the cosmological matter 
density is fixed at = 0.27. 

• The M 200 and C200 values obtained for the CFASH clus¬ 
ters using GR, Cubic Galileon and Nonlocal gravity agree 
to better than 5%, which is much smaller than the ~ 50% — 
80% precision allowed by the data at the lcr level (cf. Ta¬ 
ble [II]). In the case of the Galileon model, this is because 
the screening mechanism suppresses the modifications to 
gravity very efficiently on the scales probed by the lensing 
data, R < f? 2 oo- Oi the case of the Nonlocal gravity model, 
there are no systematic shifts in the values of M 200 and C200 


relative to those in ACDM because the boost in the gravi¬ 
tational strength is not strong enough at the redshift of the 
CFASH clusters, z ~ 0.2 — 0.9. 

• The practically unmodified lensing masses in the Galileon 
model have interesting implications for tests of gravity that 
are designed to detect differences between lensing and dy¬ 
namical mass estimates QMS1. These tests were first 
put forward in the context of models like f(R) and DGP 
that modify the dynamical potential (probed by, e.g., in¬ 
falling galaxies outside f? 2 oo)> but not the lensing potential. 
Our results show that, although the Galileon model also 
modifies the lensing potential, this does not translate into 
modified lensing masses because of the screening. How¬ 
ever, outside i? 200 > the force profile of the CFASH clus¬ 
ters in the Galileon model can be 10 — 40% higher than in 
ACDM (cf. Fig. |9j, which can affect the velocity distribu¬ 
tion of infalling galaxies, and hence, the dynamical mass 
estimates. The picture is therefore qualitatively similar to 
that of f{R) and DGP gravity, and as a result, the tech¬ 
niques of Refs. iT6]TT8l can also be applied to models like 
the Galileon model studied here. 

The existence of screening mechanisms in modified grav¬ 
ity models (such as the Galileon, which leads to practically 
unmodified lensing masses) motivates research into the lens¬ 
ing effects associated with cosmic voids. There, the density 
is low, and as a result, the fifth force effects are manifested 
more prominently. The lensing signal associated with voids 
has been detected recently in Refs. EUEa by stacking 
voids found in the galaxy distribution of Sloan Digital Sky 
Survey catalogues (see also Ref. 1123ft for an earlier fore¬ 
cast study). In the context of modified gravity. Ref. 021 
showed that voids found in simulations are on average emp¬ 
tier in f(R) gravity than in ACDM. This happens because 
the enhanced gravity facilitates the pile up of matter in the 
surrounding walls and filaments, leaving less mass inside the 
void. This translates into a stronger signature in the lensing 
signal from voids. By the same reasons, the expectation is 
that voids should also be emptier in the Cubic Galileon and 
Nonlocal gravity models. However, in these models, one has 
also the effects of the modified lensing potential, which should 
amplify the size of an eventual signature for modified gravity. 
This suggests that lensing by voids could become a very pow¬ 
erful tool to test the law of gravity outside of the Solar System. 
Such an investigation is the subject of ongoing work cm 
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Appendix A: Cluster redshifts and best-fitting parameters in 
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TABLE II. Best-fitting M 200 [W 15 Mq /hj and C200 values for all the CLASH clusters in the ACDM, Galileon and Nonlocal gravity models. 
The errors quoted correspond to the marginalized l-cr limits. 


Cluster name 
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